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ABSTRACT 

Solar Energetic Particles (SEPs) show a rich variety of spectra and relative 
abundances of many ionic species and their isotopes. A long standing puzzle 
has been the extreme enrichments of 3 He ions. The most extreme enrichments 
are observed in low fluence, the so-called impulsive, events which are believed 
to be produced at the flare site in the solar corona with little scattering and 
acceleration during transport to the Earth. In such events 3 He ions show a 
characteristic concave curved spectra in a log-log plot. In two earlier papers 
(Liu et al. 2004 and 2006) we showed how such extreme enrichments and such 
spectra can result in the model developed by Petrosian & Liu (2004), where ions 
are accelerated stochastically by plasma waves or turbulence. In this paper we 
address the relative distributions of the fluences of 3 He and 4 He ions presented 
by Ho et al. (2005) which show that while the distribution of 4 He fluence (which 
we believe is a good measure of the flare strength) like many other extensive 
characteristics of solar flare, is fairly broad, the 3 He fluence is limited to a narrow 
range. Moreover, the ratio of the fluences shows a strong correlation with the 
4 He fluence. One of the predictions of our model presented in the 2006 paper was 
presence of steep variation of the fluence ratio with the level of turbulence or the 
rate of acceleration. We show here that this feature of the model can reproduce 
the observed distribution of the fluences with very few free parameters. The 
primary reason for the success of the model in both fronts is because fully ionized 
3 He ion, with its unique charge to mass ratio, can resonantly interact with more 
plasma modes and accelerate more readily than 4 He . Essentially in most flares, 
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all background 3 He ions are accelerated to few MeV/nucleon range, while this 
happens for 4 He ions only in very strong events. A much smaller fraction of 4 He 
ions reach such energies in weaker events. 

Subject headings: Sun: flares, particle emissions - acceleration-plasmas - turbulence- 
waves 



1. INTRODUCTION 

Solar flares are excellent particle accelerators. Some of these particles on open field lines 
are observed as solar energetic particles (SEPs) at one AU or produce type III and other 
radio radiation. Those on closed field lines can be observed by the radiation they produce 
as they interact with solar plasma and fields. Electrons produce nonthermal bremsstrahlung 
and synchrotron photons in the hard X-ray and microwave range, while protons (and other 
ions) excite nuclear lines in the 1 to 7 MeV range or may produce higher energy gamma-rays 
via 7T° production and its decay. It appears that stochastic acceleration (SA) of particles by 
plasma waves or turbulence plays an important role in production of high energy particles 
and consequent plasma heating in solar flares (e.g., Ramaty 1979; Mobius et al. 1980, 1982; 
Hamilton & Petrosian 1992; Miller et al. 1997; Petrosian & Liu 2004, hereafter PL04). 
This theory was applied to the acceleration of nonthermal electrons (Miller & Ramaty 1987; 
Hamilton & Petrosian 1992). It appears that it can produce many of the observed radiative 
signatures such as broad band spectral features (Park, Petrosian & Schwartz 1997; PL04) 
and the commonly observed hard X-ray emission from the tops of flaring loops (Masuda 
et al. 1994; Petrosian & Donaghy 1999). It is also commonly believed that the observed 
relative abundances of ions in SEPs favor a SA model (e.g. Mason et al. 1986 and Mazur 
et al. 1992). More recent observations have confirmed this picture (see Mason et al. 2000, 
2002, Reames et al. 1994 and 1997, and Miller 2003). One of the most vexing problem of 
SEPs has been the enhancement of 3 He in the so-called impulsive or 3 He -rich events, which 
sometimes can be 3 — 4 orders of mag oitude above the photospherie value! There have 
been many attempts to explain this enhancement. Most of the proposed models, except the 
Ramaty and Kozlovsky (1974) model based on spalation (which has many problems), rely 
on resonant wave-particle interactions and the unique charge-to- mass ratio of 3 He (see e.g. 



addition there is charge-to-mass ratio dependent enhancement relative to the photospheric values of 
heavy ions in SEPs , and in few flares gamma-ray line emissions also points to anomalous abundance pattern 
of the accelerated ions (Share & Murphy 1998; Hua, Ramaty & Lingenfelter 1989). We will not be dealing 
with these anomalies in this paper. 
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Ibragimov & Kocharov 1977; Fisk 1978; Temerin & Roth 1992; Miller & Vinas 1993; Zhang 
1995; Paesold, Kallenbach & Benz 2003). Most of these model assume presence of some 
particular kind of waves which preferentially heats 3 He ions to a higher temperature than 
4 He ions, which then become seeds for subsequent acceleration by some (usually) unspecified 
mechanism (for more detailed discussion see Petrosian 2008). None of these earlier works 
did a compare model spectra with observations. 




Fig. 1. Left: Variation of the ratio of 3 He to 4 He fiuences with the fluence of 4 He showing a continuum of enrichments 

and a strong anti correlation. Middle: 3 He vs 4 He fiuences showing a much larger range for the latter while the former seems 
to be limited to a small range. Note that the 3 He fiuences do not concentrate at the lower end which would be the case if 
observational threshold was affecting their distribution. Right: The distribution of fiuences of 3 He and 4 He . Note that the 
high end of the 4 He distribution may be truncated because of the threshold of the fluence ratio (missing point in the lower left 
triangle of the middle panel) [From Ho05]. The fiuences are in units of particles/ (cm 2 sr MeV/nucleon). 

In two more recent papers Liu, Petrosian & Mason 2004 and 2006 (LPM04, LPM06) have 
demonstrated that a SA model by parallel propagating waves can explain both the extreme 
enhancement of 3 He and can reproduce the observed 3 He and 4 He spectra. In LPM06 
it was shown that the relative fiuences of these ions, and to a lesser extent their spectral 
indexes, depend on several model parameters so that in a large sample of events one would 
expect some dispersion in the distributions of fiuences and spectra. Ho et al. (2005; Ho05) 
analyzed a large sample of events and provide distributions of 3 He and 4 He fiuences and 
the correlations between them. Our aim here is to explore the possibility of explaining these 
observations by the above mentioned dependence of the fiuences on the model parameters. 
In particular we would like to explain the observations reproduced in Figured] which shows a 
strong anti correlation of 3 He / 4 He ratio with 4 He fluence (left panel), but shows essentially 
no correlation between the two fiuences (middle panel). More strikingly, the 3 He fluence 
distribution appears to be relatively narrow and follows a log-normal distribution, while 
4 He distribution is much broader and may have a power law distribution in the middle of 
the range, where the observational selection effects are unimportant. Often the SEPs are 
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divided into two classes; impulsive-high enrichment and gradual- normal abundance classes. 
However, as evident from the left panel of the above figure there is a continuum of enrichment 
extending over many orders of magnitude^] 

In the next section we describe some of the model characteristics that can explain these 
observations and in §3 we compare the model predictions with the observations, specifically 
the distributions of the fluences. A brief summary and conclusion is given in §4. 



2. MODEL CHARACTERISTICS 

The model used in LPM04 and LPM06 which successfully described the enrichment 
and spectra in several flares has several free parameters. As usual we have the plasma 
parameters density n, temperature T and magnetic field Bq. It turns out that the final 
results are insensitive to the temperature as long as it is higher than 2 x 10 6 K (see Fig. H] 
below), which is the case for flaring coronal loops. It also turns out that only a combination 
of density and magnetic field (y/n/B ) comes into play. We express this as the ratio of plasma 
to gyro-frequency of electrons, a = u pe /Q e which is related to the Alfven velocity in unit of 
speed of light; Pa = 5 1//2 /ct, where 5 = m e /m p is the ratio of the electron to proton masses. 
So in reality we have only one effective free plasma parameter a or Pa- On the other hand, 
several parameters are required to describe the spectrum of the turbulence. Following the 
above papers we assume broken power laws for the two relevant modes, the proton cyclotron 
(PC) and He cyclotron (HeC), with an inertial range k min < k < /c max , and similar power law 
indexes q and qh in and beyond the inertial range, respectively^ The only difference between 
the two branches is that the wave numbers fc max and k m[n for the PC mode are two times 
higher than those for the HeC mode. Finally there is the most important parameter related 
to the total energy density of turbulence, £ tot , which determines both the rate of acceleration 
and, when integrated over the volume of the source region, determines the intensity or the 
strength of the event. This parameter is the characteristic time scale r p or its inverse the 
rate defined as (see, e.g. Pryadko & Petrosian 1997) 
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with £ n - {q ~ l)£ ^ 



(^minC/^e) ^ 



9 



2 The He distribution shows a weak sign of bi modality but this is not statistically significant. In this 
paper we will ignore this feature. 

3 In LPM06 we also have an index qi describing the power law below the inertial range which is of minor 
consequence. For all practical purposes we can assume a sharp cutoff below A; m i„ which means qi — > oo. 
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for each mode. The factor of 4 arises from having two branches (PC and HeC) and two 
propagation directions of the waves (see LPM06 for details). 

As shown in LPM04 and LPM06 papers the main difference between the acceleration 
process of 3 He and 4 He is in the difference between their acceleration rate or timescales 
(r a ). The other relevant timescales, namely the loss (ti oss ) and escape (T esc ) times are essen- 
tially identical for the two ions (e.g. see left panel of Fig. 7 of LPM06). The acceleration 
timescales are different mainly at low energies (typically below one MeV/nucleon), where 
the acceleration time of 4 He is a longer (by one to two orders of magnitude). As a result at 
these low energies the 4 He acceleration time may be comparable or longer than the loss time 
which makes it difficult to accelerate 4 He ions. Most of 4 He ions are piled up below some 
energy (roughly where r a = ti oss ) and only a few of them accelerate into the observable range 
(e.g. see right panel of Fig. 7 of LPM06). However, because the acceleration times scale as 
t p while the loss time does not, for higher level of turbulence (larger £ ), the acceleration 
time may fall below the loss time so that 4 He ions can be then accelerated more readily (see 
Fig. [3] below). On the other hand, essentially independent of values of any of the above 
parameters, the 3 He acceleration time at all energies, in particular at low energies, is always 
far below its loss time so that in all cases (except for very high densities or very low values 
of T^ 1 ) 3 He ions are accelerated easily to high energies. The relative values of the escape 
and acceleration times (for both ions) determine their high energy spectral cutoffs. 

Figure [2] shows variation with energy of acceleration times of 3 He (thick lines) and 4 He 
(thin lines) and their dependence on parameters k m [ n , a and q. The remaining parameters 
qh and A; max only affect the slope of the low energy end of 4 He which does not affect the 
spectra noticeably. It is evident that the general behavior of the acceleration time scales 
described above (consisting of a low and a high energy monotonically increasing branches 
with a declining transition in between) is present in all models. These features change 
only quantitatively and often by small amounts. As expected lowering k min decreases the 
acceleration times at the high energy branch (left panel). This is because the lower k m [ n 
waves interact resonantly with higher energy ions. On the other hand, a lower value of a 
(or larger Alfven velocity or magnetization) decreases the times at the low energy branch 
(middle panel). Steeper spectra in the inertial range, produce a higher rate of acceleration 
(larger £ ; see eq. [TJ) and decrease the overall acceleration time scales (right panel) 

Note that in this and subsequent figures, fc max is in units of Vt p /c so that fc max = 2a6~ 1 / 2 
in the labels means an actual fc max = 2Q p /va = \j2j3 p /r g ^ where (3 P = 2{v t h, P /vA) 2 is 
the plasma beta, and v t h, P = ksT '/m p and r 9iP = v t h, P /Q P are the proton thermal velocity 
and gyro radius. The scale of /c max is clearly beyond the MHD regime (where the wave 
frequency to = VAk <C Q p ), but is below the proton gyro radius for the chosen parameters 
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Fig. 2. Dependence of the acceleration time of 4 He (thin, blue) and 3 He (thick, red) on k m i n (left), a (middle) and q 

(right). The lines are labeled with the corresponding numbers of each parameter. In each case the solid lines are for the fiducial 
model with a = 0.5, fc m ax = 2a<5 -1 / 2 = 10fc m i n , q = 2 and q^ = 4. 

(^/W P ~ 0.03). 

Using these acceleration rates we calculate spectra of the two ions (as in LPM04 and 
LPM06) for a range of parameters. Figure [3] shows three sets of spectra where we vary 
fcmm, a and r" 1 . In each panel the sold lines are for the fiducial model (a = 0.5, fc max = 
2a5~ 1 / 2 = I0k min ,q = 2 and = 4) chosen to fit the spectra observed by ACE/ULEIS 
for 30 Sep. 1999 event. The spectral variations here reflect the above described variations 
of the acceleration timescales. Lower fc mm (or larger inertial range) yields a larger tail for 
both ions (left panel). Variation of a has a similar and smaller effect on 3 He spectra but it 
affects the 4 He spectra dramatically; for a ~ 1 essentially there is no 4 He acceleration but 
the a ~ 1/4 model accelerates a large number of 4 He ions beyond 0.1 MeV/nucleon and 
into the observable range (middle panel). This effect is even more pronounced for increasing 
values of t~ x , where a factor of few increase in the general rate of acceleration (or the level 
of turbulence) causes a large increase of the fluence of 4 He (right panel), because, as stated 
above, its acceleration time becomes shorter than its loss time even al low energies. All these 
spectra show the same general characteristic features. While most 3 He ions are accelerated 
to high energies for essentially all model parameters appropriate for solar coronal conditions 
and reasonable level of turbulence, He ions show a characteristic lower energy bump with 
a nonthermal hard tail. In general, the lower energy bump is below the observation range 
except for low a and high values of r" 1 . Since a high level of turbulence is expected for 
brighter and stronger events, this means that we get smaller 3 He / 4 He flux or fluence ratios 
for brighter events. Note that the spectra in such cases may not agree with observations 
but this is not troublesome, because as is well established, the stronger events (the so-called 
gradual events) are associated with CMEs and shocks which most likely will modify the 
above spectra which are those of ions escaping the corona. Thus, the higher energy bumps 
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in the spectra shown here should be considered as seeds for such further acceleration during 
the transport from the lower corona to the Earth, which becomes more likely, and is expected 
to change the above spectra more significantly, for more energetic events. Thus if we give 
up the idea that there are two distinct classes of SEPs (impulsive and highly enriched and 
gradual and normal abundance) but that there is a continuum of events, which observations 
in Figure [1] show, then the above scenario implies that the main acceleration occurs in the 
solar corona. Subsequent interactions in CME shocks mainly modify the seed population 
escaping the turbulent coronal site. 




icr E io'' i io icr E icr' i 10 icr E kt 1 I 10 

E (MeV/Nucleon) E (MeV/Nucleon) E (MeV/Nucleon) 



Fie. 3. Dependence of the accelerated spectra of ^He and ^He on fc m i n (left), a. (middle) and t p 1 (right; t p 1 in units 
of r p 1,0 = 0.0055 s _1 ). The lines are labeled with the corresponding numbers of each parameter. In each case the solid lines 
are for the fiducial model with ct = 0.5, &max — 

q = 2 and q^ = 4 that is chosen to fit the data point shown 
for the 30 Sep. 1999 event observed by ACE. Note that for a better indication what energy particles dominate the spectra in a 
log-log plot we plot particle energy time fluence. 

From the spectra we can calculate the ratio of 3 He to 4 He fluences for different models 
which could be then compared with the observed ratios shown in Figure [TJ Inspection of 
observed spectra indicate that a representative ion energy would be 1 MeV/nucleon. In 
Figure H] we show the variation of this ratio with temperature (left panel) and r~ x (middle 
and right panels) for several values of other important parameters. As evident this ratio 
is most sensitive to the value of r~ l which represents the general rate of acceleration or 
the level of turbulence. The ratio can change from the highest observed value (~ 30) to 
near photospheric value (~ 2 x 10~ 4 ) for only a factor of 30 change in r" 1 . It is natural to 
expect higher level of turbulence generation (i.e. a larger value of r^ 1 ) in stronger events. 
Therefore, this predicted correlation is in agreement with the general trend of observation 
shown in Figured] (left panel), if the strength of an event is measured by the observed fluence 
of 4 He ions and most other ions like carbon, nitrogen and oxygen0 This seems reasonable 



4 It should be noted that while the observations are for fluences integrated from 0.2 to 2.0 MeV/nucleon 
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and calls for more quantitative comparison with observations and model prediction. In the 
next section we present one such comparison. 




10"' 10° 10 1 10" a 10" 1 10" a 10' 



Fig. 4. Variation of the accelerated 3 He to 4 He fluence ratio (at E = 1 MeV/nucleon) with background plasma 

temperature T (left) and Tp 1 (middle and right) for several values of other specified model parameters. The lines are labeled 
with the corresponding numbers of each parameter. In each case the open circle stands for the model that fits spectra of the 
30 Sep. 1999 event, and the solid lines are for the fiducial model with a = 0.5, fc ma x 

= 2a5~ 1 / 2 = 10fc min ,c? = 2 and q h = 4. 

Note the weak dependence on the temperature for T > 2 X 10 6 K and a strong dependence on Tp 1 for all model parameters 
with saturates at chromospheric values of the ratio. The horizontal dot-dash line shows the highest ratio observed so far (see 
Fig. Q] left). 



3. Distributions of Fluences 

We have seen that the general observed behavior of the the ratio of the fluences defined 
as R = F3/F4 is similar to the model predictions. In this section we try to put this result 
on a firmer quantitative footing by considering the observed distributions of the fluences 
of both ions as shown in Figured] (right panel). Except for the minor truncation at high 
values of F 4 , the fluence of 4 He, the observed distribution of F 3 , the fluences of 3 He, seem 
to be almost bias free and not affected significantly by the observational selection effects. 
For example, there are well defined and steep decline both at the high and low fluences 
away from the peak 3 He value of F ~ 10 3 ' 7 — , ^w C , les , — 7. This is not what one would 

J 1 u cm^sr(MeV/nucleon) 

expect if the data suffered truncation due to a low observation threshold. In such a case 
one would observe a distribution increasing up to the threshold followed by a rapid cutoff 
below it. Our model results described above also seem to predict the observed behavior. As 
stressed in previous section, the 3 He spectra and fluxes appear to be fairly independent of 
model parameters because essentially under all conditions most 3 He ions are accelerated and 



our theoretical ratios are calculated at 1 MeV/ nucleon which is near the geometric or algebraic mean of the 
range. 
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form a characteristic concave spectrum. Thus we believe that it is safe to assume that the 
observed 3 He distribution is a true representations of the intrinsic distribution (as produced 
on the Sun). This distribution can be fitted very nicely with a log-normal expression^ If we 
define the logs of the fluences and their ratio as 

LF 3 = ln(F 3 /Fo), LF 4 = ln(F 4 /F ) LR = InR, (2) 

then from fitting the observed distribution of 3 He by a log-normal form we get: 

( LF \ 2 

if) 3 (LF 3 )=<p Q exp[—^) with (T 3 ~0.22, (3) 



CT3 

which is shown on the right panel of Figure El 

Using this distribution we now derive the distribution of 4 He fluences, ip^LF^). For 
this we use the model predicted relationship between the two fluences as shown in Figure [4] 
above. We will use the two panels of this figure showing the dependence of the log of the 
fluence ratio LR on r" 1 . It turns out that most of these curves can be fitted by a simple 
function: 

HR/R Q ) = LR — lmRo = - - A — y (4) 

ln(r p 1 /r p0 ) 

The left panel of Figure [5] shows fits to the curves in the right panel of Figure H] with the 
indicated values of the the fitting parameters A, Rq and t~q (which is not the same as the 
Q = 0.0055s -1 in Fig. [3]). We shall use this relation to transfer the 3 He fluences and 
distributions to those of 4 He . 

For a given value of r' 1 the number of events with 4 He log-fluences between LF4 
and LF4 + d(LFi) (i.e. ip^LF^d^LF^)) is equal to ip 3 (LF 3 )d(LF 3 ), the number of events 
with 3 He log-fluence LF 3 = LF 4 + LR^- 1 ) and LF 3 + d(LF 3 ), whered{LF 3 ) = d(LF 4 )and 
LR{t~ 1 ) = IiiRq + A/ln(T~ l /T^). Thus we have 

/ LF 4 + LRir- 1 )^ 
MLF 4 ) = tfj 3 (LF A + LRir- 1 ]) = o exp I p M (5) 

However, we expect not a single value for r" 1 , which as stated above is a proxy for the 
strength of the event, but a broad distribution of events with different strengths, say /(t^ 1 ). 



5 The truncation shown by the shaded area in the middle panel of Figure [T] introduces a slight bias against 
detection of low fluences. We estimate that, because there are fewer events at the high 4 He fluence end, 
this means a 10 to 20% underestimation of the distribution of the three lowest values of the 3 He histogram 
(right panel, Fig. [T]). We will ignore this small correction, whose main effect is to increase the value of 03 
by a small amount. 
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Fig. 5. Left: A simple analytic fit (curves) to the model relations (points) between the fluence ratio and the acceleration 

rate or event strength represented by Tp for the three different values of fc m i n of the right panel of Figure [4] with the indicated 
fitting parameters. Right: The fitted log-normal distribution to the 3 He fluences and predicted 4 He distributions of three 
models compared with observations. The solid line which gives the best fit is for n = 2, fc m i n = 0.1fc max , the dashed line is for 
n = 2, fc m i n = 0.2fc max and the dash-dot line is for n = 1.5, fc m i n = 0.2fc max . 

Since, as argued above, the 3 He fluence distribution ^ 3 (LF 3 ) is independent of r" 1 , then for 
a population of events we have 

/do / LF + LR(t~ 1 )\ 2 

0oexp( 4 ^ P ) fir-^dr- 1 . (6) 

Every term in the above equations is determined by observations and our models except 
the distribution /(t~ , which is a reflection of the level of the distribution of the level of 
turbulence and, when multiplied by the volume of the turbulent acceleration region (which 
does not affect the 3 He / 4 He ratio), is related the overall strength of the event. Observations 
of solar flares show that most extensive characteristics which are a good measure of the 
flare strength or magnitude, such as X-ray, optical or radio fluxes, appear to obey a steep 
power law distribution, usually expressed as a cumulative distribution $(> F») oc F^ n (or 
differential distribution <f>(Fi) oc F~ n_1 ) with typically n ~ 1.5 (see, e.g. Dennis 1985 and 
reference therein). Such a distribution seems to roughly agree with the prediction of the 
so-called avalanche model proposed by Lu & Hamilton (1991). Now assuming that t~ 1 also 
obeys such a power law distribution (i.e. f(r~ l ) oc (TjT 1 )~ ( -' 1+1 - ) ) we can write the distribution 
of 4 He as: 

ip 4 (LF 4 ) = j o expl '—j e nx dx, with x = ln(r p 1 /r p0 1 ) . (7) 

Using the above relations we have calculated the 4 He fluence distribution. The results for 
three models are compared with the observations on the right panel of Figure [51 Given the 
other model parameters (k min , a etc.) we have only one free parameter namely the index n for 
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this fit. The solid line obtained for the top curve of the left panel (k min = 0.2/c max , a = 0.5), 
and for n = 2 provides a good fit to the observed distribution of 4 He fluences. In order to 
demonstrate the sensitivity of the results to the parameters we also show two other model 
predictions based on slightly different parameter values. These results provide additional 
quantitative evidence (beside those given in LPM04 and LPM06) on the validity of the SA 
of SEPs by turbulence, and indicate that with this kind of analysis one can begin to constrain 
model parameters. 



4. SUMMARY AND CONCLUSIONS 

In this paper we have carried out further comparison between the prediction of models 
based on stochastic acceleration of SEP ions by turbulence. In our earlier works (LPM04, 
LPM06) we demonstrated that the extreme enrichments of 3 He and spectra of 4 He and 3 He 
observed in several events can be naturally described in such a model. Using the results based 
on this model, here we consider the relative distributions of 4 He and 3 He fluences derived 
from a large sample of event by Ho05. We show that with some simple and reasonable 
assumptions we can explain the general features of these observations as well. 

These are clearly preliminary results and are intended to demonstrate that in addition 
to modeling only few bright events it is also important to look at population as a whole and 
ascertain that a model which can explain the detail characteristics of individual events can 
also agree with the distributions of observables for a large sample of events. Here we have 
shown how the dispersion in one parameter, namely the acceleration rate or the strength of 
the flare, can account for the observed distributions of fluences. The key assumption here 
is that the amount of produced turbulence (represented by T p v ) has a wide dispersion and 
obeys a power law distribution similar to that observed for other extensive parameters that 
give a measure of the strength of a flare. The dispersion in other model parameters can 
also influence the final outcome. However, the dispersion of most of the other important 
parameters, like intensive parameters temperature, density and magnetic field, are expected 
to be much smaller than that of an extensive parameter like the overall strength of the 
event, the amount of turbulence produced, the flare volume etc. In addition, as shown in the 
previous section, the intensive parameters play a lesser role than the extensive parameter 
Tp 1 in determining the relative characteristics of 3 He and 4 He . Given the dispersion of 
any other parameter one may carry out similar integration over its range. However, for the 
reasons given above we expect smaller changes in the shapes of predicted distribution due 
to dispersion of most of the intensive parameters. Given a more extensive set of data such 
improvements may be needed can be carried out. 
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The existing data may be used to test some of our assumptions, in particular the as- 
sumption of constancy of the 3 He distribution. We intend to address these in future works. 
We can also make the above results more robust by using model fluences integrated over 
the same spectral range as the observations instead of fluences at 1 MeV/nucleon. One can 
also expand this approach and address the distributions of other characteristics besides the 
fluence, sue as the spectral indexes or break energies (if any). The available data contain 
this information but require more analysis. 
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